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Abstract 

We study BPS saturated objects with axial geometry (wall junctions, vor- 
tices) in generalized Wess-Zumino models. It is observed that the tension of 
such objects is negative in general (although "exceptional" models are pos- 
sible). We show how an ambiguity in the definition of central charges does 
not affect physical quantities, and we comment on the stability of the junc- 
tions and vortices. We illustrate these issues in two classes of models with Zn 
symmetry. On the basis of analytical large N calculations and numerical cal- 
culations at finite N, we argue that the domain wall junctions in these models 
are indeed BPS saturated, and we calculate the junction tensions explicitly. 
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1 Introduction 



Studying the BPS domain walls in supersymmetric theories is interesting, especially 
in gauge theories Q, because one can get information on non-trivial dynamical 
features. For example, it was found that the tension of such domain walls determines 
the decay rate of "false "vacua in large N supersymmetric gluodynamics @, |[]. 

It was noted 0] that theories with either a U(l) or global symmetry may 
contain BPS objects with axial geometry. Vortex- like objects which conserve 1/4 of 
the original supersymmetry were considered in |J, while domain wall junctions of 
the hub and spoke type, which also preserve 1/4 of supersymmetry, were discussed 
in |Q, |5], [7| . General considerations of the tensorial central charges responsible for the 
saturation of the vortices and junctions were presented in [Q, |9], [10| . The existence of 
BPS wall junctions also leads to non-trivial dynamical consequences. For instance, 
in [|ll|] the large N behavior of the BPS junction tension was used to argue in 
favor of soliton-like heavy hadrons with mass proportional to iV in supersymmetric 
gluodynamics. 

The study of possible dynamical implementation of the BPS wall junctions has 
just begun. Apart from the original works in which the existence of such junctions 



was first noted, they were discussed in [|iT| , [12], [13], [14| . In [|TT] a generalized Wess- 



Zumino model, which we will refer to as Model I, with one chiral superfield X and 
the superpotential 

was considered. This model is related to the Veneziano-Yankielowicz []L5[] effective 
Lagrangian, which, in turn, models supersymmetric Yang-Mills theory. The param- 
eter N is an integer that corresponds to the the number of colors of the gauge group 
SU(iV) of SUSY gluodynamics. The superpotential (1) was originally suggested for 
that purpose in |16|| . At any given N it possesses a Z^ symmetry which is spon- 



taneously broken. In the large N limit an analytic solution for the BPS domain 



walls was found in fPT |. The most essential feature of the solution is the fact that 
the wall width scales as iV -1 , and vanishes at iV — > oo. This result, in conjunction 
with two BPS junction configurations presented in at N = oo, implies that BPS 
junctions do exist at large N in the theory with superpotential (1). The junction 
tension was found to scale as N 2 at large N. 



In |L3[] the first analytic solution for a BPS junction was found in a specific 
generalized Wess-Zumino model. Among stimulating findings in this work is the 
fact that the junction tension turned out to be negative in this model. The model 
has Z 3 symmetry. It is derived from a SU(2) Yang-Mills theory with extended 
supersymmetry (A/ - = 2) and one matter flavor perturbed by an adjoint scalar mass. 
The original model contains three pairs of chiral superfields and, in addition, one 
extra chiral superfield. In fact, for the purpose of studying the BPS walls and 
wall junctions, the model can be further simplified. It can be easily stripped of 
its inessential features and cast into the form of a Z^ symmetric model. (N is an 
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arbitrary integer, not necessarily equal to 3). The generalized and simplified model, 
which we refer to as Model II, contains N chiral superfields M k with k = 1,2, ...,N 
and another chiral superfield T. The superpotential is 



W = N<-A 



+ BT\, (2) 



£ (T - Nq^e^) Ml 
Lfe=l 

where q = e l2lT / N . The constants A and B are assumed to be N independent and 
positive. 

Model I is invariant under the transformation 

X(x, 6) -> e 2mk/N X(x, 6') , & = e - mk/N 9 ; (3) 

the corresponding transformation of the superpotential is 

W(X(x, 9)) -> e 2mk/N W{X{x, 6')) . (4) 

Model II is invariant under the transformation 

T(x, 6) -> e 2 " fc/JV T(x, 0') , Mt(x, 6) -> M,_ fe (x, 0') 0' = e-" 1 ^^ ; (5) 

the corresponding transformation of the superpotential is 

W(T(x, 9), M(x, 6)) -> e 2mfc/7V W(X(x, 0'), 0')) • (6) 

Model I is among the simplest supersymmetric models allowing one to study 
BPS junctions and calculate tensions. [] The model has iV distinct vacua in which 
the scalar component of X takes the vacuum expectation value j\f e t7T ( 2k - 3 )/ N ) with 
k — 1, 2, .., iV. The perturbative spectrum consists of a chiral multiplet with mass 
m = N. 

Let us now briefly discuss the vacuum structure and spectrum of Model II. In 
this model there are 2N distinct, physically equivalent vacua. In these vacua the 
scalar components of the T and M fields take the vacuum expectation values 

T = Ne in{2k - 3)/N , M k = ±J^, M, = for l^k, (7) 

V 

with k — 1,2, .., N. The perturbative spectrum contains two chiral multiplets, linear 
combinations of the field T and the field that obtains a vacuum expectation 
value, with mass m 12 = 2\J ABN, and (N — 1) chiral multiplets with mass 

m k = AN 2 A(l-cosTr(k-l)/N), (8) 

for k = 2,3,..,N. 



x It is worth noticing that in a different context its 2-dimensional reduction was studied in |18[ 
where it was shown to be integrable. 
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Apart from the symmetry, Eq.(|5]), Model II has additional Z^ symmetry, 
since any of the fields can be transformed as M\. — > —M^. The superpotential 
in Eq.(^) is obviously invariant under this transformation. In our consideration of 
the wall junctions below we will disregard the Z2 symmetry, limiting ourselves just 
to a Zn family of vacua. (This corresponds to the choice of the + sign in front of 



A/B in Eq.(|7|)). The walls and junctions of the type where (some of) the vacua 
lie outside this family may or may not be BPS saturated. We leave this issue open 
for future investigations. 

Needless to say that in both models I and II the walls connecting the distinct 
vacua and the corresponding junctions are topologically stable. (The issue of sta- 
bility of the BPS junctions is discussed separately in Sect.(^).) They need not be 
BPS saturated, generally speaking. We know, empirically, that sometimes in certain 
models equations of the BPS saturation have no solutions with the proper boundary 



conditions |18], [L9| even in the case of isolated walls, let alone the wall junctions. 

It was known previously that the isolated walls in Model I are saturated, see 
|fT6| . In Appendix A we give an analytical solution to the BPS equation for the 
domain walls in Model II for a specific ratio of the parameters A and B\ we have 
verified numerically that these walls also saturate the BPS bounds for other values 
of the parameters. In Sect.(|6|) we provide evidence that the domain wall junctions 
in Models I and II are also BPS saturated. Were the walls and junctions non- 
saturated, one could still calculate the relevant central charges. They would then 
represent lower bounds on the corresponding tensions. 

Most considerations in this paper concern the basic domain wall junction, the 
junction with iV sectors in a model with Zn symmetry. But in Appendix B we 
present an analytic solution to the BPS equation for the class of triple junctions in 
Model II for a specific ratio of the parameters A and B, and iV a multiple of three. 
For N = 3, the basic junction and the triple junction coincide. Non-basic junctions 



in Model I are studied in fL4 |. Networks of domain walls in Model I are studied in 



In this paper we consider BPS wall junctions and calculate their tensions in 
Models I and II specified above. We also consider general aspects, including the 
sign of the tension, the ambiguity of the central charges, and stability, that are valid 
beyond Model I and II. These models are used to illustrate the general assertions. 

The organization of our paper is as follows. In Sect.(|2|) we describe a general 
method to calculate the BPS bound on the tension of domain wall junctions. We 
show that the junction tension is typically negative, although we do not exclude 
exotic models with positive tension. We then use this method to explicitly calculate 
the junction tension in Models I and II. (Sees. 3 and 4). For Model I we calculate 
the junction tension analytically in next to leading order in 1/N, and numerically 
for finite values of N. For Model II the junction tension is calculated analytically for 
all values of N and all values of the parameters A and B. In Sect.(|5|) we show that 
ambiguities in the central charges cancel in the calculation of the junction tension, 
and we illustrate this cancellation by recalculating the tension in Model II with a 
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different definition of the central charges. In Sect.(|!]) we demonstrate by numerical 
- and in certain limits, analytical - calculations that in both models I and II the 
wall junctions are, in fact, saturated. In Sect.(|7|) we comment on the stability of the 
BPS junctions. The appendices contain analytical results for domain walls and the 
triple junction in Model II. 



2 Generalities of the solutions with the axial ge- 
ometry 

For both Models I and II it is possible to calculate, using the general theory of 
the central charges, the BPS bounds on the junction tension. The general formula 
worked out in QTOI implies 



M 

T 



- j> a k dx k + 2 I dn k S k , (9) 



where a k is an axial current. At the classical level only the component of the current 
built from the scalar fields is important. Moreover, 

Si, 2 = {ReS, ImS} , (10) 

where £ is related to the superpotential. The integrals run over the large circle (see 
Fig. (|l])) with R — > oo, and dx (dn) is an infinitessimal vector along (perpendicular 
to) the contour. The precise definitions of a k and E have (correlated) ambiguities. 



The sum in Eq.(|9]) is unambiguous, however [[10]. The above ambiguity can be 
exploited to optimize the calculation of the junction tension. This will be explained 
below in two particular applications, Model I and II. 

We will consider the wall junction of the "hub and spoke" type. The following 
conventions are convenient (although not necessary). We will orient the wall spokes 
as indicated in Fig. ([I]), namely, the hub is at the origin, the "first" spoke runs along 
the x axis in the positive direction, the "second" runs at an angle 2-k/N ', and so on. 
In the point P the theory "sits" in the first vacuum, in the point Q in the second, 
etc. This configuration is topologically stable. 

We will assume that the Zn symmetry is realized through multiplication of (some 
of) the fields by a phase. This is certainly true in Models I and II. If on the junction 
solution (the lowest components of) the fields X, T (generically, $) rotate as 

we will say that the rotation weight of the field $ is equal to r. Configurations with 
the positive rotation weight will be referred to as junctions (vortices), while those 
with the negative weight will be referred to as anti-junctions (anti- vortices). It is 
worth reminding that 



( = x + iy, d c = -(d x -id y ). (12) 
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Figure 1: A "standard" convention regarding the wall junction. The large circle of 
radius R enters the definition of E(R) while the circle with the segments removed is 
relevant to E(R). The box of size L x L around the point G is blown up in Fig.(Q). 
L is an auxiliary parameter chosen to satisfy the constraint / <C L <C R. 
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Figure 2: A blow up of the box in Fig.(|l|). The shaded rectangle represents the wall, 
the length of the interval GG' is of order L 2 /R. 
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The BPS wall junctions satisfy the equations 



29 <*' = e aI7- (13) 



where the phase 5 depends on the phase of the superpotential (see below). The 
equations for ant i-j unctions is 



i.e. ( is replaced by (. Thus, in our convention the phase of $ grows in the anticlock- 
wise directions for the junction, and in the clock-wise direction for the anti-junction. 

The superpotential can always be chosen in such a way that in the vacuum P 
its phase is exp(— iir/N) while in the vacuum Q it is exp(iir/N). See Fig.(|l|) for the 
definition of P and Q. (Note that Re{W} is then positive in both vacua). Then the 
phase 5 in Eq. ( |T3|) must be set to zero in both Models I and II. 

For the wall junctions the rotation weights of X and T are obviously equal to 1. 

Now, we proceed to discuss general features of the tension associated with wall 
junctions. In Fig. ([I]) the energy of the junction configuration (per unit length) is 
defined as the integral over the volume energy density over the area inside the circle, 
where it is assumed that the circle radius R tends to infinity, 

M r 

E(R) = - r = H(x,y)dxdy = T 1 R + T 2 + ..., R -> oo . (15) 

length J\f\<R 

Here the dots denote terms vanishing in the limit R — > oo. It is tempting to say that 
T 2 is associated with the (1/2,1/2) central charge while T\ with the (1,0) central 
charge (see Eq.(|), the first and second term, respectively). We hasten to note that 
these central charges, being considered individually, are ambiguous [10]. It is only 



the combination in Eq.(g) which is fixed unambiguously. It is intuitively clear that 
T\ = iVT wall where T wall is the tension of the isolated wall. Moreover, it is also clear 
that T 2 (which will be also referred to as Taction) is typically negative. Indeed, from 
Eq.fllSD it follows that, for small R, E(R) = cR 2 , where c is a numerical coefficient. 
At large R the quadratic dependence on R changes to linear, see Eq. (p~5f) . Matching 
these two scaling laws straightforwardly, connecting the parabola and the straight 
line at the value of R where their slope is equal, we conclude that T 2 is forced to 
be negative. Models where it could be positive would require, in essence, two scales 
and an intermediate regime between the parabolic and linear regime. In order to 
illustrate the typical situation, we plot the energy E as a function of R for Model 
II with N = 3, B = V3 and A = B/9 in Fig.®. 

Actually, this figure gives the energy E defined as the integral over T~C(x, y) over 
the equilateral triangle with distance R from the center to each of the sides instead of 
the energy in a circle. The orientation of the triangle is such that each of the sides 
is perpendicular to a spoke. We made this modification because for this triangle 



6 



E 




Figure 3: Energy E as a function of radius R (dotted line) for the triple junction of 
model II with N = 3, A = B/N 2 and B = \/3, compared to the line E = T X R + T 2 
(solid line). This example illustrates the typical junction with negative tension. 



the energy converges to the linear relation Eq.(|15D exponentially fast, whereas for 
a circle the departure from the linear behavior is of the order 1/R. (As explained 
below, the coefficient of the 1/R term is just a geometrical factor taking into account 
the curvature of the circle as it intersects the domain wall.) 

For the particular values of the parameters that we chose, there is an analytical 
wall junction solution of the BPS equations, which is presented in Appendix B. 
The energy E(R) can therefore be calculated exactly and is shown in Fig.(Q). It is 
clear from Fig.(|^) that even though there is a change from a convex to a concave 
regime, E(R) converges to the linear regime so fast that Tj unct i on is negative. The 
tensions can be calculated analytically and take the values T\ = 162 (see Sect.(|6|)) 
and T 2 = — 27v^3/2 (see Sect.(|J)). Similar graphs for Model I junctions are presented 
in Ref.||14||. For comparison, we sketch what a similar graph would look like for an 



exotic model with positive junction tension in Fig.(f|). 

Below we give a general prove of the fact that Ti = iVT wall , and illustrate more 
quantitatively that Tj unction is negative in both Models I and II, and in other "natu- 
ral" models. Of crucial importance is the fact that the wall width (i.e. the transverse 
dimension inside which the energy density is nonvanishing, while outside it vanishes 
with exponential accuracy) is finite (i.e. R independent). This width is denoted 
by I, see Figs. (P and (|). In both Models I and II, I is actually 0(l/N). At the 
moment we do not assume iV to be large. It can be any integer > 2. If R — > oo, one 
can replace the integration domain in Eq. flT5l) - the interior of the circle can be re- 
placed by the interior of the circle with the segments removed. If E(R) is defined as 
the integral over the interior of the circle with the segments removed, the difference 
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Figure 4: Sketch of the energy E as a function of radius R for a hypothetical, exotic 
model with the positive junction tension. 



E(R) — E(R) tends to zero in the limit R — > oo. (The length of the straight line 
DH is L, where L is an arbitrary parameter, / <C L <C R.) This can be clearly seen 
from Fig.(0), which presents a detail of the domain near the point G. The area of 
the removed segment is 0(L 3 /R). The parameter L does not scale with R when R 
tends to oo, hence the area of the removed segment tends to zero. Since the volume 
energy density is finite, the difference E(R) — E(R) ~ R" 1 . 

For the BPS saturated configuration the energy E(R) or E(R) can be rewritten 
as a combination of the corresponding contour integrals. The master equation for 
the wall junction is presented in Eq.([|). It is assumed in this equation that the phase 
5 in Eq. ( fI3[ ) is zero, which is the case with our convention. There is an ambiguity 
in S and a k in this master formula: a certain reshuffling is possible due to the 
ambiguity in the definition of the supercurrent and the energy-momentum tensor. 
This is explained in great detail in Ref. |10||. We will return to the discussion of this 



ambiguity later on. For the time being let us stick to the canonic definitions. Then 

E = W , S ii2 = {ReW, ImW} , (16) 



and 



\ ]T ($*z3*^ - QiidkQi) = i £ (W S h) , (17) 



where k — 1,2. The sum runs over all scalar fields involved in the solution. Fields 
which are purely real obviously drop out. Only complex fields contribute; for such 
fields, if $ = pe irj , 

a k = p 2 d k rj. (18) 
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Note that Eq. (|ITD for the current appearing in the master equation is valid for the 
wall junctions. For the anti-junctions one must reverse the sign in Eqs. (|T7|) and 



d y 3> = i— . (19) 



How do we see that with our convention the phase 5 in Eq. ( |13|) is zero? To this 
end, let us consider this equation far away from the "hub" , on the line DH. On this 
line Eq. (|T3|) becomes 

. dW 
~W 

Since AW = W(Q) — W(P) = 2i\W\ sin(7r/A r ) in our convention, this is exactly 
the BPS equation for the (isolated) wall, as it should be since the equation for the 
junction should go into the equation for the wall far away from the "hub." 

To begin with, we will calculate the tension T\. Since due to the Zn symmetry, 
the contributions of each of the N sectors in Fig. ([I]) are identical, it is sufficient to 
calculate the integral corresponding to one sector, 



T X R = 2N / dn k S k , (20) 

JPDHQ 

where on PD and HQ the integral runs along the circle while on DH it runs along 
the straight line. On PD and HQ the superpotential W assumes the vacuum values 

W = |W| vac e™ /7V on HQ , W = | WU c e~™ /N on PD , 

and does not change (with exponential accuracy). Therefore, 

/ dn k S k = R\W\ vac If 1 cos ( ^ + a) + I*'" cos ( ^ - a) \ , (21) 

JPD+HQ \J-tt/N \N ) J-y \N J 



where the angle 7 is defined in Fig.(^). On the straight line DH we have 

-L/2 

/DH J-L/2 



/ dn k S k = / dyReW. (22) 

JDH J-L/2 



With the exponential accuracy the fields on the straight line DH satisfy the BPS 
equations for the isolated wall, see (0). The solution is known to have an integral 
of motion (see e.g. 0), which with our phase convention is nothing but Re{W}. 
Thus, Re{W} does not change along the line DH and reduces to its vacuum value, 
|VV| vac cos(7r/iV), which is one and the same in both vacua, P and Q. As a result, 
the sum of two integrals (|2l|) and (^) is 



/ dn k S k = it!|WU c 2sin - 7 )+L |W| vac cos ^ = 2i?|W| vac sin ^-+0(R- 1 ) 

JPDHQ \JM / JM JM 

(23) 



'PDHQ 

Moreover, 



2|W| vac sin^ = \W{P)-W{Q)\ = ^T wall , (24) 
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where T wal i stands for the tension of the isolated BPS wall. Combining Eqs. (f20|), 
(p3|) and fl2"4|) we arrive at 

Ti = iVT wall , (25) 

q.e.d. 

Now we pass to the discussion of the tension T 2 . According to Eq.([|) for the 
wall junction tension, T 2 is given by the integral over the large circle over 

T 2 = - <j> dx k a k , (26) 

where a k is given in Eq. ( |T7j ) or Eq. ([Tj^). Again, we take into account the fact 
that all N sectors give one and the same contribution, and we will integrate over 
one sector only, starting from P and ending at Q. Then 

T 2 = -N [ dx k a k . (27) 

JPDHQ 

On the intervals PD and HQ, the fields assume their vacuum values. They do 
not change (with exponential accuracy), and therefore, a k = 0, so that there is no 
contribution to T 2 . A nonvanishing contribution comes from the DH segment of the 
integral. On this segment we can disregard the "hub" and other "spokes" of the 
junction, considering the relevant part of the wall as that of the isolated BPS wall. 

In Model I, there is only one field, X. In Model II there are N + 1 fields, T 
and Mi. However, for the BPS wall, M/s assume real values (see Appendix A) and, 
thus, the fields Mg do not contribute. Let us parametrize the fields X and T as pe lv . 
In both cases rj(P) = —tt/N and i](Q) = tt/N. Then 

T 2 = -N [ dyp 2 (x,y)d y r ] (x,y) = -N(p 2 ) ( V (Q) - V (P)) = -2tt( P 2 ) . (28) 

J DH 

Here (p 2 ) is the average value of the modulus of the field on the solution on the 
segment DH, i.e. it is the average over the wall. First, we notice that since (p 2 ) 
is positive, the value of the junction tension Tj unction is negative. The only way to 
escape this, is to have some (complex) fields involved in the junction which rotate in 
the unnatural direction (i.e. while the "natural" fields in the junction rotate in the 
anticlock-wise direction, this "exotic" field must rotate clock-wise). Rotation means 
here winding of the phase, from smaller to larger values. In addition, the "exotic" 
field must have the average value of the modulus larger than that of the natural 
fields, so that it can overcome the negative contribution to Tj unct i on coming from the 
"natural" fields. In all models considered so far of which we are aware, the "exotic" 
fields do not appear .0 

Second, we note that (p 2 ) inside the wall does not differ too much from its 
vacuum value. For instance, in the Model II (see Appendix A) the vacuum value 

2 Let us note in passing that similar consideration determining the sign of the tension T 2 is 
applicable to the U{\) vortices of the type considered in S. 



10 



1^1 vac = N 2 while in the middle of the wall |T| middle = N 2 [cos(ir / N)] 2 . The average 
value (p 2 ) lies between these two extremes. At large N the average value tends to 
its vacuum one, (p 2 ) = N 2 + 0(1). 

Theorem For BPS junctions Tj unct i n is determined, through the master equa- 
tion, by the solution for the isolated wall. Unlike T\, the domain wall junction is 
non-holomorphic, and, generally speaking, depends on the details of the wall solu- 
tion. In the limit N —>■ oo the formula for Ij unc tion becomes universal. 

Proof 

The first part of the theorem is proved by the consideration preceding Eq. (E^) . 
The last part of the theorem will be established below (see also |TT| ). 



3 Taction from the isolated wall in Model I 

First, we will illustrate how this works in Model I. Let us consider this model in 



the limit of large N. In this limit the wall junction tension was calculated in flTOf . 
^junction = —2irN 2 . The corresponding calculation is trivial. We are now interested 
in the leading 1/N correction. It can be found using the solution for the isolated 
BPS wall obtained in [17]. The field X is parametrized as 

X = N(l-£) e^ N , (29) 

where a(y) is a function of y which falls off from <r* in the middle of the wall (y = 0) 
to zero at y = oo. Here ex* is the positive root of the equation 

(1 - o*)e a * = -1 , (7*« 1.278. (30) 

The phase r changes from zero in the middle of the wall {y = 0) to tt at y = oo, 

t(ct) = 7i — arccos ((1 — a) exp a) , (31) 



where we modified the solution presented in |L7] to bring it in accordance with our 
convention. 

Now in the segment under consideration 

dr 

a y = (N- 2ct) — + 0{1/N) , (32) 
dy 

and 

- / a v dy = -2tcN + 4 / dya— = -2vriV + 4 / r(a)da 
jdh ' Jo dy Jo 

= -2vriv(l-^ + 0(l/iV 2 )) . (33) 
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This implies, in turn that 



T junction = -2ixN 2 (l - ^ + 0(l/iV 2 )) . (34) 

The tension in Model I can also be calculated numerically for any value of N. In 
Fig.(^) we compare the numerically determined tension with the leading and next 
to leading analytic approximations for N = 2 k and k = 2, 10. The value of N 
therefore ranges between 4 and 1024. A numerical profile of the domain wall was 
obtained by simulation of the second order equation of motion for the field X on a 
discrete lattice, using a forward predicting algorithm. The lattice spacing was taken 
much smaller than the width of the wall, and at the same time the size of the lattice 
was chosen much larger than this width. The value of the field was fixed to take on 
the two vacuum values at the two ends of the one dimensional lattice. A dissipation 
term was added to the equation of motion so that the field relaxes to the minimum 
energy configuration. An interpolating profile between the two vacuum expectation 
values was chosen as the initial condition. Errors due to the finite lattice spacing 
are well under control because the shape of the wall profile is smooth. Errors due 
to the finite size of the lattice are suppressed because the wall profile tends to the 
vacuum expectation values exponentially fast. After the field had come to rest, the 
junction tension was then obtained from Eq.fl2"T|) by numerical integration. 

It may be worth a few words to motivate our choice of the second order equation 
of motion over the first order BPS equation to obtain the numerical domain wall 
profile. Simulation of the equation of motion yields a domain wall profile even if the 
wall is not BPS saturated, but in the case of Model I this does not offer an advantage, 
as the walls are known to be BPS saturated. However, numerical integration of the 
first order BPS equation is fraught with an instability. This can best be seen from 
an approximate solution of the BPS equation near a vacuum value of the field. 
Parametrizing the field as X = N(e m ^ N + a + bi), the solution to the linearized BPS 
equation for a and b (valid for \a\, \b\ -C 1) is 

71 -N 2 ii 71 N 2 ii 

a = c\ cos — e y + c 2 cos — e y 

b = Cl (l- s in^)e-^-c 2 (l + sin^)e^, (35) 

where c\ and c 2 are real constants. For a wall connecting the vacuum with k — 1 at 
y = — oo to the vacuum with k — 2 at y — oo, c 2 = 0. The value of c\ is associated 
with the location of the center of the wall. In simple numerical procedures to 
integrate the BPS equation, c 2 will be small, but not equal to zero. As a result, the 
numerical solution will show runaway behavior far away from the center of the wall. 



12 



T / N 

2 




10 k, N=2 



Figure 5: The tension of the junction in Model I as a function of N. The dashed 
line reflects the leading order contribution in 1/N. The solid line includes the next 
to leading order contribution, and the dots represent the numerical calculation. 



4 Taction in Model II 

Although in general Tj unction depends on the detailed shape of the domain walls 
surrounding the junction, this is not the case for Model II. In fact, the tension of 
the junction in Model II does not depend on A and B and can be calculated exactly 
even though an analytic solution for the domain wall is not available for generic 
ratios of A and B. The only assumption is that the fields Mi and M 2 are real on 
the domain wall profile. This is consistent with the explicit domain wall solution 
for the special ratio of A and B discussed in Appendix A, and it is confirmed by 
numerical solutions of the domain wall profile for generic ratios. The field T then 
takes the form T = —it + N cosir/N on the domain wall, where t is a real function 
with asymptotic values t(— oo) = iVsin7r/iV and t(oo) = —Nsinn/N. The junction 
tension is now equal to 

//•OO (j^ 
dx k a k = N N costt/N— dx = -N 3 sm2ir/N. (36) 
J-oo dx 

The tension therefore only depends on the asymptotic values of the field T on the 
domain wall and is independent of the detailed shape of the wall. In this sense, the 
junction tension Tj unction in Model II is similar to the (1, 0) tension of domain walls. 
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5 Reshuffling between the (1,0) and (1/2, 1/2) cen- 
tral charges and the unique expression for T\ 2 



As was discussed at length in [10], in the case of wall junctions in supersymmetric 
theories, one deals, in fact, with two central charges - one appears in the anticom- 
mutator {QuQp} (this is the (1, 0) central charge), another appears in the anticom- 
mutator {Q a Qp\ (this is the (1/2, 1/2) central charge). The individual expressions 
for these central charges are not unique. They depend on the expression for the 
supercurrent one starts with. This is due to the fact that the effect under discussion 
is subtle - we are interested in integrals over total derivatives. The supercurrent has 
an ambiguity. Say, in the Wess-Zumino model with one chiral superfield one can 
add to the supercurrent the term 



V2 



9 a p (^ f ) + d p$ tyatf) - 3ep a dj (^ 7 0t)] , (37) 



which is conserved (nondynamically) and presents a full derivative. It has no impact 
on the supercharges defined as 



Q a =fd"xJ° a , J£ = I (an $ " J aP p (3* 



2 



(for definitions see ||). The contribution of AJ^ in Q a is of the form / d 3 xV(...). 
Adding the term in Eq. fl37|) one changes both central charges. 

If we do not add the term AJ Q/3/ j in the supercurrent at all (by the way, this 
will correspond to the canonic energy- momentum tensor), then we get Eq.(BI) with 



S and a defined in Eqs. (|1J) and fll7[) . Assume that we add to the supercurrent the 
term AJ a ^ as in Eq. fl3"7p, with the coefficient indicated in this expression. What 
changes? 

It is not difficult to show |j| that the changes are as follows: (ii) the energy- 
momentum tensor appearing in the anticommutator of two supercharges acquires a 
full- derivative (conserved) term built from the scalar fields; this term is such that 
9^ now vanishes for purely cubic superpotential, i.e. for the conformally invariant 
theory; (ii) £ in Eq. ( |16|) now becomes 

1 dW 

X = W--$— , S li2 = {ReE,Im£}; (39) 

(iii) simultaneously, the expression — § a^dx^ in Eq.(||) changes too. Let us keep, for 
simplicity, the definition of the current ( |T7| ) intact. Then a "new" master equation 
takes the form 

M If f 

length = ~3 J °' k<i ' Xk + 2 J dnkSk ' ( 40 ) 

where S is defined in Eq 
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It is instructive to check that, in spite of the reshuffling in the (1,0) and (1/2, 1/2) 
central charges, the final results for the tensions T 12 remain unambiguous. They are 
insensitive to this reshuffling. One can do this check explicitly using the analytic 
solution for the BPS wall which exist in Model II, see Appendix. (A general proof 
was presented in [fL0| .) As previously, we split the contour integral in N sectors; 
we will explicitly consider only the transition PQ (Fig. ([]])). On the segments PD 
and HQ the extra term in the new master equation J2^dW/d^ vanishes, since 
these segments correspond to vacua, where <9W/<9$ = 0. It is only the interval DH 
that contributes to the integral over (—2/3) <&dW/d<&. One can readily convince 
oneself that the contribution due to M 12 appears only in ImS = and, thus, 
drops out. Only the field T is relevant. Due to the fact that only the segments of 
the type DH contribute, it is immediately clear that the above extra term in the 
new master equation changes only the tension T 2 leaving T\ intact. 

Now, at R — > oo one can replace the DH segment by that for an isolated BPS 
wall, using explicit expressions collected in Appendix A. In this way we get 



2 /" L<9Wl 2 Ar2 . 



2tt 
~N 



(41) 



and, consequently, 



2 / dn k AS k = -^N 3 sin ^ = 5 2 (T 2 ) , (42) 



where 



3 <9$ y ' v 3 <9$ 
Using the same approach we readily find that the term 



j> dx k a k 



-N 3 sin ^ ; (44) 



this term was in the old master equation. In the new master equation we have one 
third of it, i.e. the contribution to T 2 coming from the axial current is one third of 
that in Eq. d§, 

5i(T 2 ) = -Vsin^. (45) 
Combining this with S 2 (T 2 ) from Eq. (|42|), we arrive at 

T 2 = 5 X (T 2 ) + 5 2 (T 2 ) = -N 3 sin ^ , (46) 

exactly the same result for Tj unc tion that was obtained from the old master equation 
where only the axial current contributed to it. The reshuffling in the central charges 
does take place; there is no impact on the physical quantities, however. 
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6 BPS saturation of junctions in Models I and II 



In this section we consider the question whether the domain wall junctions in Models 
I and II saturate the BPS bound. For certain limiting values of the parameters, 
the BPS equations for domain wall junctions in Models I and II allow an analytic 



solution. This is the case for the N — ► oo limit in model I, as was shown in |L0 
And in Appendix B we present an analytic solution for the triple junction in Model 
II, which can be obtained when N is a multiple of 3 and for a particular ratio of A 
and B. This junction is a generalization of the analytic solution in jnj. For N = 3 
the triple junction is also the basic junction. The above values of the parameters 
seem only special because they allow an analytical solution, and not for any other 
reason. Moreover, for any value of the parameters there are sensible approximate 
solutions to the BPS equation far away from the center of the junction (domain 
walls surrounding the center) and near the center (string-like solutions). All this 
information strongly suggests that the junctions are BPS saturated for all values of 
the parameters. 

In order to further test this hypothesis, we performed a numerical analysis. We 
simulated the second order equations of motion for both Models I and II on a lattice 
using a simple forward predicting algorithm. The procedure is a two-dimensional 
generalization of the procedure described in Sect.(|3]) to generate numerical domain 
wall profiles. The lattice spacing was chosen to be much smaller than the charac- 
teristic size I of the junctions, and at the same time the size of the lattice was much 
larger than I. We put the appropriate domain wall profiles (numerically generated) 
on the edges of the lattice as boundary conditions. A dissipation term was added to 
the equations of motion, so that the field configuration relaxes to the configuration 
with lowest energy consistent with the boundary conditions. 

We then determined whether the numerically obtained junction profile saturated 
the BPS bound by comparing the energy of the configuration to E(R) = T\R + T 2 . 
For Model I, T\ = 4A 4 /(iV + 1) shi7r/iV, and T 2 was determined numerically as in 
Sect. (|). For Model II, T\ = AN 3 B sin n/N and T 2 is given in Eq.(||). Numerical 
errors due to the finite lattice spacing were well under control because the junction 
shape is smooth, and errors due to the finite size of the lattice were suppressed 
because the junction profile tends toward the domain wall solutions exponentially 
fast away from its center. 

We performed this procedure for Model I with N = 4, and for Model II with 
N = A and various values of the parameters A and B. In each case we found that the 
junction was indeed BPS saturated. This fact, in conjunction with the analytical 
results, lends strong support to the conjecture that the basic junctions we have 
discussed in Models I and II are BPS saturated for all values of the parameters. 
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7 Comments on Stability of the BPS wall junc- 
tions 



A general question which deserves a brief discussion is the stability (local and 
global) of the BPS wall junctions. In this issue two features are crucial: (i) the 
two-dimensional manifold perpendicular to the junction axis (the plane) is noncom- 
pact, i.e. the junction "spokes" extend to infinity; (ii) the wall junction is not a 
two-dimensional soliton; rather, each "spoke" represents a wall which extends in the 
direction perpendicular to the plane. This extension is infinite. 

As was discussed above, the mass of the BPS wall junction can be expressed in 
terms of two central charges which reduce, in turn, to the contour integral over the 
large circle (see Fig.([|a), where we consider a Z% configuration, for definiteness) . 
This means that any localized perturbation of the junction configuration, both in 
the plane (Fig||b), and in the perpendicular direction, will lead to an encreased 
energy, i.e. the BPS wall junction is stable against such localized perturbations. 

A global perturbation, when two "spokes" approach each other and then even- 
tually fuse into one, is energetically favorable (Figs. (|5].c) and (E|.d). This is due to 
the fact that the isolated BPS wall tension for the transition from the vacuum I to 
the vacuum J has the form 

T 13 = 2|Wj - Wi\ = 2|W vac | x \e 2 ™ J ' N - e ^ iI/N \ , (47) 

in the models under consideration, and, hence 

Ti_ n + 7n_m > 7i_ m . (48) 

In other words, eating up an intermediate vacuum we lower the energy of the con- 
figuration. 

Were the wall sizes finite, the configuration with six spokes in Fig.(|.a) could 
quantum-mechanically tunnel into the configuration with five spokes (Fig.(|].d), un- 
der the barrier. However, since the walls are infinite in both dimensions, this cannot 
happen. The boundary conditions at infinities make the BPS wall junctions stable 
both classically and quantum-mechanically. For further consideration see Ref. |L4 



8 Conclusions 

In supersymmetric models with the spontaneously broken Zn symmetry we devel- 
oped techniques for (i) establishing the BPS nature of the wall junctions of the "hub 
and spokes" type and for (ii) calculating the tensions 7\ and Tj unction . It is shown 
that in "natural" models Tj unction < 0. The conditions under which the "unnatural" 
sign Tj unc tion > can occur are discussed. 
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Figure 6: Stability of the wall junction, a.) Regular BPS junction in a model with 
Zq. b.) A localized deformation of the wall junction, c.) Squeezing out the vacuum 
II. d.) The wall junction with the vacuum II missing. 
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Appendix A 

Here we present an analytic solution for the BPS wall in Model II which can be 
readily found for a specific ratio of the coefficients A and B, namely 

A ~ 2- (49) 



[2iVsin(7r/iV)] 

For what follows it is convenient to introduce the parameter A which determines the 
width of the wall, 

A = r 1 = - , f — . (50) 
2sin(7r/iV) v 1 

Note that at large N the parameter A scales as N; correspondingly, / ~ N^ 1 . 
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We will consider the wall interpolating between the vacuum 1 ( "initial" vacuum) 
at y — — oo and the vacuum 2 ("final" vacuum) at y — +00. 
The vacuum 1 lies at 

IB 

T = Ne~ wi/N , M 1 = J—, M 2 = M 3 = .... = 0. (51) 

V 

The vacuum 2 lies at 

T = Ne wi/N , M 2 = , Mi = M 3 = .... = . (52) 

For the superpotential in the vacua 1 and 2 we have 

W(l) = BN 2 e~™ /N , W(2) = BN 2 e wi/N , (53) 

so that 

AW = 2iBN 2 sm(n/N) . (54) 
Correspondingly, the BPS equations for the fields T, Mi j2 have the form 

T = t^l , M E = i~ . (55) 



Here the dot denotes differentiation over y. The solution of these equations is 

7T / . 7T 

iv +2iV l sin iv 



T = A^cos|- + ^fsin|-)tanh(Ay), (56) 



and 



/S 1 -tanh(Ay) „, B 1 + tanh(Ay) 
Ml = VI 2~^' M2 = 2~^' (57) 



Appendix B 

Here we present an analytic solution for BPS walls and triple junctions in Model 
II, with N a multiple of 3. An analytic solution for the triple junction exists if the 
coefficients A and B satisfy the relation 

A = m- m 

Note that the ratio of A and B depends on N, in contrast with what was assumed 
in the main text. For this particular ratio of A and B, the parameter A which 
determines the width of the wall and the size of the junction takes the form 
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We will consider a junction that interpolates between three of the N vacua. Vacuum 
1 lies at 

T = Ne' m ' N ) M 1 = ^|, M k = 0,k^l, (60) 

vacuum 2 lies at 

T = Ne 2m/3 e - m /N ^ Mjy/3+i = ^| ) Mk = 0,k^N/3 + l, (61) 

and vacuum 3 lies at 

T = iVe 47 ™/ 3 e-™ /7V , M 2N/3+1 = ^| , M k = , k ^ 2N/3 + 1. (62) 

Let us first focus on a domain wall interpolating between vacua 1 and 2 along the y 
axis. Such a wall is a solution to the BPS equations 

T = ie V- M '= fe "ajar- M "'— "saw < 63 » 

where the dot indicates differentiation with respect to y. When the phase takes the 
value 5 = 7r/3 — n/N, the wall configuration tends to vacuum 1 at y — > — oo and to 
vacuum 2 at y — > oo. The analytic solution to the BPS equation is 

T = ^ e ~ m/N (e™ /3 - v / 3e- m / 6 tanh(Ay)) , (64) 



and 

[B 1 - tanh(Ay) /j? 1 + tanh(Ay) 
1 = ^A 2 ' M7V / 3+1 = V 1 * ' ^ ' 

Having the solution for the wall in hand, we now continue with the triple junction. 
In accordance with the convention in the main text, our choice of 5 implies that 
the spoke in between the first and the second vacuum points in the direction of the 
positive x axis. For N = 3, the phase 5 = 0, as in the main text. For other values 
of N, insisting on 5 = would lead to a very awkward orientation of the junction. 
The junction is the solution to the BPS equations 

and 

n dM x _ udW n dM N/3+1 _ lS dW dM 2N/3+1 _^ iS dW 



d( 9Mi' ' d( dM N/3+1 ' ' d( dM 2N/3+1 

(67) 



20 



with ( = x + iy. Denning 

g( x , y) = e -My-x/Vz) + e Ky+*/Vz) + e -2Xx/V~^ (6g) 

the triple junction BPS configuration takes the form 

-X(y-x/V3) I p 2TTi/3 p \{y+x/V3) i p Ani/ 3 p -2\x/V3 

T = Ne- m / N - — , — , (69) 

9{x,y) 

and 

e -\(y-x/V3) fg e \(y+x/V3) fg e ~2\x/V3 

' ' V A g(x,y) ' N/3+1 V A g(x,y) ' 27V/3+1 V A g(x,y) 

. (70) 

For x — > oo the solution tends to the domain wall configuration along the y axis 
which was discussed before. For x — > — oo and y = the field configuration tends 
towards vacuum 3. The tension for the triple junction is 

T 2 = ~V3N 2 . (71) 

This is an exact result, valid for any value of A and B, not just the special ratio 
for which the analytic solution for the junction is obtained. For N = 3 this result 
agrees with the tension presented in Eq . fl46|) in the main text. 
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